Hexagonal ABC semiconductors as ferroelectrics 
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We use a first-principles rational- design ap- 
proach to identify a previously-unrecognized class 
of ferroelectric materials in the PQ^mc LiGaGe 
structure type. We calculate structural param- 
eters, polarization and ferroelectric well depths 
both for reported and as-yet hypothetical repre- 
sentatives of this class. Our results provide guid- 
ance for the experimental realization and further 
investigation of high-performance materials suit- 
able for practical applications. 

A rapidly developing paradigm for the rational design 
of functional materials is based on the first-principles 
study of large families of known and as yet unreported 
compounds. First-principles calculations of structure and 
properties are used first to explore the microscopic origins 
and establish design principles for the functional proper- 
ties of interest, and then to screen a large number of both 
equilibrium and metastable phases to identify promising 
candidate systems [TR3]. One recent study showed the 
semiconducting members of the ABC half-Heusler fam- 
ily to be piezoelectric, with a range of piezoelectric prop- 
erties comparable to the much-studied ^4i303 perovskite 
oxides [4]. 

A ferroelectric is a material with a polar phase pro- 
duced by a structural transition from a nonpolar high- 
symmetry paraelectric state, with an electric polarization 
that can be switched between two or more symmetry- 
related variants by application of an electric field [5]. 
The rational design of new ferroelectric materials is moti- 
vated both by fundamental scientific interest and by po- 
tential technological applications [6]. New materials can 
offer better performance, including reduction in switch- 
ing time, in coercive field and in fatigue, operation at 
higher or lower temperatures, and the possibility of bet- 
ter integration with other materials based on structural 
or chemical compatibility. New ferroelectrics with lower 
band gaps for photoactive applications [7H§] are also of 
interest. Additional practical advantages could include 
decreased toxicity, for example Pb-free [ID] , and possible 
multifunctionality. 

Any polar structure (if insulating) could potentially 
support ferroelectricity if the barrier to switching is low 
enough [TTH13] . We therefore can search for new ferro- 
electric semiconductors by targeting intermetallic com- 
pounds in polar space groups and screening both reported 
and hypothetical compounds to find insulating represen- 
tatives with a low barrier to uniform switching through 



a nonpolar reference phase, which provides an indication 
of the barrier to realistic switching. ABC compounds 
with polar space group P&^mc in the LiGaGe structure 
type [T1HT7] are a promising target class. This structure, 
shown in Figure [T] is a hexagonal variant of the half- 
Heusler structure and can be described as a wurtzite 
structure "stuffed" with a third cation [18]. The Inor- 
ganic Crystal Structural Database (ICSD) |19j includes 
18 ABC compounds in this structure type that do not 
contain an /-block element. We can classify these com- 
binations into the following groups: I-III-IV (LiGaGe), 
I-II-V (LiBeSb), I-XII-V (LiZnSb), XI-III-IV (CuYSn), 
XI-II-V (AgCaBi) and II-XII-IV (CaZnSn). In addition, 
we find that six entries (CuScSn, CuYSn, AuYSn, Ag- 
CaBi, CaZnSn and CaHgSn) are also reported with non- 
polar Pdz/mmc symmetry in the ZrBeSi structure type, 
which we identify as the nonpolar reference phase. It has 
been previously noted that the P&s/mmc ZrBeSi struc- 
ture can be obtained by a symmetry-restoring distortion 
of the LiGaGe structure in which the buckling of the 
atomic planes in the wurtzite structure is eliminated; this 
relationship is analogous to that of the wurtzite structure 
to the metastable hexagonal structure of ScN [20J [21] . 

In this paper, we use first-principles methods to estab- 
lish a new class of ferroelectrics in the LiGaGe structure 
type and to identify promising candidate materials for 
further investigation. Specifically, we compute the struc- 
tural parameters, band gap, polarization, and barrier to 
uniform switching of the eighteen reported and 70 as- 
yet-hypothetical ABC compounds in the LiGaGe struc- 
ture type. We identify several insulating combinations 
with polarization comparable to or greater than that of 
BaTiOs, and uniform switching barriers comparable to 
or less than that of PbTi03. For all insulating combi- 
nations studied, we find that the band gaps are in the 
semiconducting range; the lower band gaps could be use- 
ful for photoactive applications [7H5]. These candidate 
ferroelectrics offer promise for experimental investigation 
and for the future development of new high-performance 
materials for practical applications. 

First principles computations were performed with the 
ABINIT package [52] . The local density approximation 
(LDA) and a 4 x 4 x 4 Monkhorst-Pack sampling of the 
Brillouin zone [53] were used for all calculations, except 
for the Berry phase polarization [24] [25] calculations, for 
which a 8 x 8 x 8 grid was used. All atoms were repre- 
sented by norm-conserving optimized [26 designed non- 
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local [57] pseudopotentials, generated with the OPIUM 
code [28]. All calculations were performed with a plane 
wave cutoff of 50 Ry. 

We first consider the eighteen non-rare-earth com- 
pounds that have been experimentally reported in the 
LiGaGe structure type in the ICSD [3T]. For each com- 
bination, we optimized the structural parameters for each 
of the three structural variants ABC, ABC, and ABC_, 
where the underscore indicates the element at Wyckoff 
position 2a, which stuffs the wurtzite structure comprised 
by the other two elements. The results for the lowest en- 
ergy structural variant are reported in Table [I] The com- 
puted structural parameters generally show good agree- 
ment with experimental values, with the underestimate 
of lattice constants characteristic of LDA calculations, 
about 1-2 % for a and as large as 3-4 % for c. 

There are only two insulating compounds in the set 
of reported LiGaGe-type compounds: LiBeSb (I-II-V) 
and LiZnSb (I-XII-V); in each the stuffing atom is the 
monovalent element Li. The total of 8 s and p valence 
electrons is expected to improve the likelihood of band 
gap formation 29 . The band gap for LiBeSb is 1.71 
eV (indirect) and for LiZnSb it is 0.67 eV (direct). The 
computed polarizations of 0.59 C/m 2 (LiBeSb) and 0.56 
C/m 2 (LiZnSb) are larger than that of BaTi03. To as- 
sess switchability, we compute the energy difference be- 
tween the polar state and the nonpolar high-symmetry 
reference state. While we recognize that ferroelectrics do 
not switch by uniform change of the polarization through 
the high-symmetry state, the energy barrier for uniform 
switching can be used to assess the possibility of realistic 
switching by comparing with the values for known fer- 
roelectrics: 0.2 eV for PbTi0 3 and 0.02 eV for BaTi0 3 
[3"D] . In the present case, this comparison suggests that 
the nominal barriers in LiBeSb and LiZnSb (0.58 eV and 
0.80 eV) are too high for switchability to be likely. 

To search for candidate LiGaGe-type ferroelectrics 
with lower barriers, we consider equiatomic combina- 
tions of three distinct constituent elements ABC with 
valences given by I-II-V or I-XII-V, with I=(Li, Na, K), 
II=(Be, Mg, Ca, Sr, Ba), XII=(Zn) and V=(P, As, Sb, 
Bi) . This generates a total of 72 candidate combinations 
to be searched, only two of which are included in Ta- 
ble I. We optimize the structural parameters for each of 
the three variants corresponding to the three choices of 
element for the 2a position. 

We find that 6 of the 72 combinations are found to 
be metallic in the lowest-energy structural variant: the 
computed structural parameters and AE, the energy rel- 
ative to the relaxed high-symmetry P63 /mmc phase, for 
these are given in Table [n] For the remaining 66 of the 
72 combinations, the predicted lowest-energy variant is 
insulating, with computed band gap ranging from 0.04 
eV to 1.81 eV; since DFT tends to underestimate band 
gaps, we expect that the actual fraction of insulating 
compounds will be slightly higher than our calculations 



would indicate. Of these 66 compounds, 49 have relaxed 
to the higher non-polar Pd^/mmc symmetry; results for 
these compounds are given in the supplementary mate- 
rial. For the 17 polar insulating compounds, we also 
compute the spontaneous polarization; results for these 



compounds are given in Table III 

Thus, we have narrowed the search for new LiGaGe- 
type ferroelectrics to seventeen polar insulating combina- 
tions. For this set, we see in Figure[2]that AE has a pos- 
itive correlation with polarization, as would be expected 
in a simple double-well model. The eight compounds 
LiBeP, LiCaBi, NaMgP, NaMgAs, NaZnSb, NaMgBi, 
KMgSb, and KMgBi have polarizations comparable to 
that of BaTi0 3 and AE <0.25 eV, in the range favor- 
able for ferroelectric switching, and therefore are promis- 
ing candidates for ferroelectricity. 

The ferroelectric double well for NaMgP, shown in the 
inset of Figure [2] is representative of this group. The 
key to the switchability of the polarization is that the 
wurtzite substructure is not characterized by ideal sp 3 
bonding like in ZnO, which would require breaking and 
reforming of rigid bonds to switch. Rather, the structure 
should be understood as the buckling of the flat planes 
of the PQz/mmc structure, with sp 2 bonding. We define 
a buckling parameter, d, as the distance along the c-axis 
between the inequivalent atoms in the buckled plane. We 
find that d decreases as the size of the stuffing ion in- 
creases from Li to Na to K, weakening the interplanar 
bond so that the barrier to switching is reduced most for 
compounds containing K. This structural trend directly 
affects the polarization, which arises from a combination 
of the buckling and the displacement of the planes rela- 
tive to the stuffing cation. Since changing the sense of 
the buckling does not involve breaking and reforming of 
bonds, the barrier to switching can be low enough for 
ferroelectricity. 

A key question is that of the prospects for synthe- 
sis of the candidate compounds in the desired struc- 
ture. Of the eight compounds we have identified as 
candidate ferroelectrics, six have reported structures in 
ICSD. Five are reported in space group PA/nmm (LiBeP, 
NaMgAs, NaZnSb, KMgSb, KMgBi), and one in space 
group Prima (LiCaBi). Of the two for which there is 
no reported structure, results from a recent theoretical 
study [3] predict them to be of PA/nmm (NaMgP) or 
P2i/c (NaMgBi) symmetry. 02 However, it could still 
be possible to synthesize at least some of our candidate 
LiGaGe-type ferroelectrics as metastable phases, in cases 
in which the LiGaGe structure type is sufficiently close 
in energy to the ground state. In particular, for NaZnSb 
the LiGaGe phase is only 0.04 eV per f.u. higher in 
energy than the P4/nmm ground state, which makes 
the metastable phase quite accessible. Furthermore, for 
NaZnSb the energy difference between the lowest energy 
variant and the next (with Zn as the stuffing atom) is 
0.51 eV, suggesting that it will be possible to obtain full 
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chemical ordering. 

If these compounds are grown as epitaxial hims, this 
would provide an additional route to engineering the po- 
larization, the switching barrier, and the relative stability 
of the LiGaGe phase. For example, first-principles cal- 
culations show that 3% tensile strain in the (0001) plane 
reduces the AE of the reported compounds LiZnP by 
0.11 cV, LiZnAs by 0.12 eV and LiMgAs by 0.08 eV, to 
0.57 eV, 0.60 eV and 0.36 eV, respectively. Strain could 
also promote a polar instability in the insulating nonpolar 
P63 / 'mmc compounds. Of the forty-nine compounds we 
have identified as nonpolar insulators, six are reported in 
ICSD with PQs/mmc symmetry (see supplemental Table 
1), and the previously mentioned theoretical study found 
five additional compounds with this structure. First prin- 
ciples calculations of the zone-center phonon frequencies 
for six selected compounds (LiBaSb, NaBeSb, NaCaBi, 
KZnAs, KZnSb, and KBaSb) show that in each case the 
frequency of the lowest frequency polar mode is below 
100 cm -1 . However, the coupling of this mode to (0001) 
epitaxial strain is not strong enough to produce an insta- 
bility in the range ± 4 % in any of the six compounds we 
tested. 

In conclusion, we have used first-principles methods to 
establish a new class of ferroelectrics in the LiGaGe struc- 
ture type and to identify promising candidate materials 
for further investigation. Through targeted synthesis, 
LiGaGe-type compounds could potentially be developed 
as a valuable class of ferroelectric and piezoelectric ma- 
terials; other structure types with substructures related 
to wurtzite could similarly yield systems with switchable 
polarization. This is a specific application of a larger- 
scale strategy to identify new ferroelectrics by targeting 
polar insulating compounds not previously recognized as 
ferroelectric and tuning the composition and other con- 
trol parameters, such as epitaxial strain, and/or modify- 
ing the structure by intercalation of atoms to reduce the 
barrier to polarization switching. The identification of 
ferroelectricity in classes of materials in which it was pre- 
viously unrecognized offers the possibility of optimizing 
properties and combining polarization with other func- 
tional properties, including magnetism, to produce multi- 
functional behavior of fundamental scientific interest and 
for groundbreaking technological applications. 
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FIG. 1: The PGs/mmc ZrBeSi structure shown in a) is 
the nonpolar high-symmetry reference structure for the polar 
PG^mc LiGaGe structure shown in b). The two structures are 
related by a buckling of the planes formed by atoms at Wyck- 
off positions 2b (dark blue) and 26' (gold) and displacements 
of the planes relative to the stuffing atom at 2a (green). 
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FIG. 2: Difference in energy between the low (PG^mc) and 
high (PQi/mmc) symmetry ABC structures vs polarization 
for all polar insulating combinations. Combinations with 
A=IA (black circles) are less likely to be switchable than those 
with yl=Na (red circles) or A=K (green circles). The two re- 
ported compounds LiBeSb and LiZnSb are labeled. The in- 
set shows the characteristic ferroelectric double well energy of 
NaMgP as a function of polar distortion obtained by linear in- 
terpolation between the polar PG^mc and nonpolar PGzjmmc 
structures. 
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a 

(A) 

4.14 
4.35 
4.48 
4.51 
4.63 



LiGaGe 
ScCuSn 
YCuSn 
YCuPb 
YAgSn 
ScAuGe 4.26 
ScAuSn 4.51 
YAuSi 
YAuGe 
YAuSn 
LiBeSb 
LiZnSb 
LiZnBi 
CaAgBi 4.73 
CaZnSn 4.60 
CaHgSn 4.75 
SrHgSn 4.88 
SrHgPb 4.89 



4.25 
4.36 
4.61 
4.09 
4.38 
4.46 



4.18 
4.39 
4.54 
4.56 
4.68 
4.31 
4.59 
4.29 
4.41 
4.64 
4.15 
4.43 
4.58 
4.81 
4.66 
4.80 
4.89 
5.00 



c 

(A) 

6.71 
6.80 
7.15 
7.19 
7.24 
6.68 
7.07 
7.36 
7.08 
7.29 
6.64 
7.08 
7.23 
7.56 
7.33 
7.55 
7.80 
8.05 



6.78 
6.83 
7.27 
7.33 
7.37 
6.85 
7.20 
7.55 
7.31 
7.37 
6.74 
7.16 
7.38 
7.83 
7.63 
7.76 
8.22 
8.17 



Z2b 

0.31 
0.33 
0.32 
0.32 
0.31 
0.30 
0.35 
0.27 
0.28 
0.33 
0.35 
0.29 
0.28 
0.32 
0.31 
0.31 
0.29 
0.29 



(0.31 
(0.33 
(0.32 
(0.32 
(0.31 
(0.30 
(0.34 
(0.27 
(0.28 
(0.32 
(0.34 
(0.30 
(0.28 
(0.31 
(0.30 
(0.30 
(0.30 
(0.30 



z 2b 

0.70 
0.73 
0.73 
0.73 
0.70 
0.71 
0.73 
0.72 
0.71 
0.73 
0.73 
0.67 
0.66 
0.72 
0.72 
0.71 
0.69 
0.69 



(0.70 
(0.73 
(0.73 
(0.73 
(0.72 
(0.70 
(0.73 
(0.73 
(0.72 
(0.73 
(0.73 
(0.69 
(0.66 
(0.72 
(0.72 
(0.72 
(0.72 
(0.72 



AE 
(eV) 
0.80 
0.45 
0.32 
0.34 
0.39 
0.61 
0.75 
0.05 
0.26 
0.70 
0.58 
0.80 
0.62 
0.36 
0.32 
0.53 
0.36 
0.19 



TABLE I: First-principles results for the eighteen LiGaGe- 
structure compounds reported in ICSD. The computed 
ground state structural parameters are compared to exper- 
imental values, given in parentheses (in some cases, a trans- 
formation has been applied to facilitate comparison). The 
origin is chosen so that the atom at position 2a (underlined) 
is at 2=0. AE is the energy difference between the polar state 
and the relaxed nonpolar high-symmetry reference state. 



ABC 

LiBeP 

LiMgP 

LiZnP 

LiBeAs 

LiMgAs 

LiZnAs 

LiBeSb 

LiZnSb 

LiBeBi 

LiCaBi 

NaMgP 

NaMgAs 

NaMgSb 

NaZnSb 

NaMgBi 

KMgSb 

KMgBi 



a 

(A) 

3.634 
4.134 
3.924 
4.091 
4.292 
4.108 
4.094 
4.376 
4.179 
4.679 
4.424 
4.549 
4.868 
4.558 
4.957 
5.017 
5.092 



c 

(A) 

5.833 
7.211 
6.365 
6.636 
7.483 
6.673 
6.645 
7.081 
6.806 
7.503 
6.877 
7.262 
7.584 
7.332 
7.584 
7.789 
8.005 



Z2b z 2b ' 



0.295 
0.413 
0.341 
0.265 
0.414 
0.335 
0.267 
0.288 
0.262 
0.287 
0.310 
0.314 
0.307 
0.285 
0.302 
0.290 
0.291 



0.686 
0.773 
0.722 
0.647 
0.774 
0.715 
0.650 
0.669 
0.643 
0.731 
0.716 
0.715 
0.705 
0.677 
0.704 
0.697 
0.698 



0.11 
0.14 
0.12 
0.12 



d £/gap 

(eV) 
1.51(i 
1.81(i 
1.33(i 
1.71(i 
0.14 1.48(d 
0.12 0.97(d 
0.12 1.13(1 
0.12 0.67(d 
0.12 0.45(i 
0.06 0.16(i 
0.09 1.17(d 
0.10 0.67(d 
0.10 0.69(d 
0.11 0.20(d 
0.10 0.14(i 
0.09 0.59(d 
0.09 0.15(i 



AE P 

(eV) (C/m 2 ) 
0.24 0.85 



0.42 
0.65 
0.30 
0.44 
0.72 
0.58 
0.80 
0.62 
0.01 
0.20 
0.23 
0.29 
0.16 
0.25 
0.08 
0.15 



0.91 
0.84 
0.67 
0.59 
0.75 
0.59 
0.56 
0.55 
0.19 
0.49 
0.48 
0.43 
0.49 
0.42 
0.40 
0.31 



TABLE III: First-principles results for the lowest-energy vari- 
ant of each of the 17 polar insulating ABC combinations in 
the LiGaGe structure type. The structure is specified by the 
lattice constants and the internal structural parameters for 
two of the atoms; the origin is chosen so that the atom at 
position 2a (underlined) is at 2=0. The buckling parameter 
d is described in the text. The computed values for the band 
gap -Eg a p (d=direct and i=indirect), AE (as in Tables I and 
II), and polarization P are included. 
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LiMgBi 


4.602 


7.051 
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0.718 


0.08 


0.42 


LiZnBi 


4.461 


7.233 


0.277 


0.658 


0.12 


0.62 


NaZnBi 


4.653 


7.479 


0.283 


0.675 


0.11 


0.13 


KSrBi 


5.757 


7.044 


0.250 


0.750 








KZnBi 


4.582 


10.135 


0.250 


0.750 









TABLE II: First-principles results for the lowest-energy vari- 
ant of each of the 6 metallic ABC combinations in the LiGaGe 
structure type. The buckling parameter d is described in the 
text. AE, as in Table 1. 
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